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We study the low energy behavior of QCD Green functions in the limit that the baryon chemical 
potential is much larger than the QCD scale parameter Aqcd- We show that there is a systematic 
low energy expansion in powers of where w is the energy and m is the screening scale. 

This expansion is valid even if the effective quark-gluon coupling g is not small. The expansion is 
purely perturbative in the magnetic regime |fc| 3>fco. If the external momenta and energies satisfy 
\k\ ~feo, planar, abelian ladder diagrams involving the full quark propagator have to be resummed 
but the corresponding Dyson-Schwinger equations are closed. 


Landau’s theory of Fermi liquids has been remarkably 
successful in many areas of physics P|. Landau argued 
that although strongly interacting many-particle systems 
are very complicated in general, the situation greatly 
simplifies if the temperature is very low. In this limit 
the properties of the system can be understood in terms 
of weakly interacting excitations called “quasi-particles”. 
The complicated dynamics of the underlying many-body 
systems can be summarized in terms of just a few pa¬ 
rameters that characterize the quasi-particles and their 
scattering amplitudes. 

In dense quark matter the presence of unscreened mag¬ 
netic interactions invalidates the basic assumptions of 
Landau Fermi liquid theory illili Weak coupling 
calculations suggest that non-Fermi liquid corrections to 
the specific heat in the normal phase ill or the gap 
in the superfluid phase @ are substantial. However, 
weak coupling methods are restricted to densities that 
are much larger than the densities that can be achieved 
deep inside a neutron star. 

This raises the question whether there is a system¬ 
atic low-energy expansion that can be used in much the 
same way that Landau theory is used in nuclear and con¬ 
densed matter physics. There has been a lot of progress 
in this direction, mostly based on effective theories of 
dense QCD l^. In this paper we show that QCD 
Green functions have low energy expansions in powers of 
where uj is the characteristic energy scale of 
the process and m is the screening scale. This expansion 
is reliable as long as w <§; m even if the gauge coupling g 
is of order one. We provide a brief exposition of our main 
results. A detailed diag rammatic analysis will appear in 
a companion paper |Il| . 

Our starting point is an effective field theory (EFT) 
for fermions near the Fermi surface interacting with un¬ 
screened gauge fields. We follow Hong and expand 
the fermionic part of the effective action in derivatives of 
the fields, 

Cf = (1) 
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Here, if±v describes particles and holes with momenta 
p= and = (1, in) labels the local Fermi veloc¬ 


ity. We will write I = lo + l\\+l± where Z|| and l± are the 
components of I parallel and orthogonal to v. We con¬ 
centrate on two patches on opposite sides of the Fermi 
surface as these are the only channels that lead to kine¬ 
matic enhancements. The effect of fluctuations above the 
cutoff A is encoded in low energy constants. The leading 
coefficients are the shift Sg. of the Fermi energy, the Fermi 
velocity vp = |n|, the renormalization factors Z^^ and Z±, 
the gauge coupling g at the scale A, and the four-fermion 
couplings and in the forward (Zero Sound) and 
back-to-back (BCS) channels. 

The dominant effect in the gluon sector is due to hard 
dense quark loops (HDLs) in gluonic correlation func¬ 
tions. These fluctuations involve small external energies 
Lo but hard loop momenta p ~ /r in a narrow interval of 
width Lo. Hard dense loops are non-analytic in the gluon 
energy and momentum and cannot be represented as lo¬ 
cal operators. However, as long as hard fluctuations do 
not change the symmetries of the groundstate all other 
corrections are analytic. We can then write the gluonic 
lagrangian as Cg = Co + Crdl + • ■. where Cq is the 
free gluon la CTang ian and £hdl is the HDL generating 
functional 
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This term describes screening and damping of soft gluon 
modes due to particle-hole pairs on the entire Fermi sur¬ 
face. In perturbation theory the dynamical gluon mass 
is given by m? = Nfg'^g ^where g is the QCD 
coupling constant. 

We now analyze the dynamics of the theory governed 
by equ. itna . We first note that particle-hole loops have 
already been integrated out and are represented by the 
HDL term. The effective theory describes the interaction 
of particles and holes with soft gluons which do not sig¬ 
nificantly change their velocity v. Since electric fields are 
screened the interaction at low energies is dominated by 
the exchange of magnetic gluons. Magnetic gluons are 
weakly damped in the kinematic regime |A:o| <C |fc|. In 
this regime 
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and we observe that the propagator becomes large if 

|fc| - > [/cqI . (4) 


This implies that the gluon is very far off its energy shell 
and not a propagating state. We will compute a general 
diagram by picking up the pole in the quark propagator, 
and integrate over the cut in the gluon propagator using 
the kinematics dictated by equ. 0- In order for a quark 
to absorb the large momentum carried by a gluon and 
stay close to the Fermi surface this momentum has to be 
transverse to the momentum of the quark. This means 
that the term fc^/(2/i) in the quark propagator is relevant 
and has to be kept at leading order. Equation 0 shows 
that (2^) fco as Aq ^0. This means that the pole of 
the quark propagator is governed by the condition A|| ^ 
A^/(2^). We conclude that quark and gluon momenta 
scale with respect to an external energy scale uj as 

Aq ^ a;, A|| ~ //J,, k±^m^uj3. (5) 


We will refer to the regime in which all momenta, includ¬ 
ing external ones, satisfy the scaling relation © as the 
magnetic regime. A similar regime was identified in the 
context of gauge theories of condensed matter systems 
la The scaling relations © are the basis of the low 
energy expansion in ultradegenerate matter. 

In the low energy regime propagators and vertices can 
be simplified even further. The quark and gluon propa¬ 
gators are 
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FIG. 1: One-loop contributions to the qnark self energy and 
the quark-gluon vertex. In the magnetic regime the graphs 
scale as tjlog(tj), and respectively. 


the expansion involves positive powers of and the 
magnetic regime is completely perturbative. 

As a simple example we consider the fermion self en¬ 
ergy. The one-loop diagram is 
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with Cf = {N^—1)/{2Nc). This expression shows a num¬ 
ber of interesting features. First we observe that the 
longitudinal and transverse momentum integrations fac¬ 
torize. The longitudinal momentum integral is performed 
by picking up the pole in the quark propagator. The re¬ 
sult is independent of the external momenta and only 
depends on the external energy. The transverse momen¬ 
tum integral is logarithmically divergent. We find ii 




and the quark gluon vertex is gZ^\Vi{X°‘/2). Higher or¬ 
der corrections can be found by expanding the quark and 
gluon propagators as well as the HDL vertices in powers 
of the small parameter e = w/m m We observe that 
the transverse projector in the gluon propagator simpli¬ 
fies because k± k\\. We also note that in the mag¬ 
netic regime the factor po in the quark propagator can 
be dropped since po^P||. 

Using these expressions we can show that the power of 
e associated with a Feynman diagram always increases 
with the number of loops and the number of higher- 
order vertices. One way to see this is to rescale the fields 
in the effective lagrangian so that the kinetic terms are 
scale invariant under the transformation 
{e~^XQ,e~'^/^x\\,e~'^/^x_i_)- The scaling behavior of the 
fields is Ip ^ e^^^ip and Ai —> We find now 

that the scaling dimension of all interaction terms is pos¬ 
itive. The quark gluon vertex scales as the HDL 
three gluon vertex scales as and both the quark- 

two-gluon vertex induced by the D\ term in equ. 
and the four gluon vertex scale as e. Since higher order 
diagrams involve at least one pair of quark gluon vertices 


where uj=po. We have absorbed the logarithmic cutoff 
dependence into the low energy constant Zy. In gen¬ 
eral this result depends on two unknown parameters g 
and As where Ae = Aexp(97r^(Z|| — I)/(5 ^Z||Uf)) and 
A = 2h.\/{nm?) is related to the transverse momentum 
cutoff. If the coupling is small, the scale is determined by 
the exchange of electric gluons and we find Ae = 2^^^™/?!. 
We observe that the self energy correction is large and 
has to be included in the propagator whenever its en¬ 
ergy dependence is relevant. We showed previously that 
rainbow diagrams do not give corrections of the form 
log(a;)” 0- Equ. © shows that hi gher order cor¬ 
rections are suppressed by powers of [l^ . 

The scaling arguments apply to arbitrary Green func¬ 
tions in the magnetic regime, see Eig.^ Exceptions occur 
if the external fields have small momenta of the order of 
the external energy scale. This situation can occur in 
quark-quark scattering amplitudes or in vertex functions 
for external currents like the weak interaction |l(i| . Con¬ 
sider the one-loop vertex correction for a color singlet 
vertex F^ = eZ\\v^. In the magnetic regime the graph 
scales like . This is confirmed by an explicit calcu- 
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FIG. 2: One-loop correction to the vertex of an external cnr- 
rent and the BCS interaction. Both diagrams are kinemati¬ 
cally enhanced and scale as log(a;) and log^(ai), respectively. 
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TABLE I: Leading order expressions for the low energy con¬ 
stants appearing in the effective Lagrangian for magnetic 
modes equ. Q in the weak coupling limit. 


lation 0. In the time-like regime po — qo > \p~^ 

r'‘(p.,) = es^c^v>t I 1 

r dk\\ „ 

X J Z^ii Sij {k + p) S-g {k + q) , (10) 


where p, q are the momenta of the external quarks. The 
important point is that if we combine the fermionic prop¬ 
agators using Fermi’s trick in order to resolve the pole in 
the longitudinal momentum integration, the large com¬ 
ponents A:|| and k‘^/{2fi) of the propagators cancel and 
the result becomes sensitive to the small scales p, q 


Sij {k+p) Sij {k+q) 


Sij{k+p) - Sff{k + q) 
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where constant factors are suppressed. As a consequence 
the result is enhanced by a factor 1/e^/^. This enhance¬ 
ment is analogous to the one occurring in the HDL case. 
In the limit po — qo ^0 the /c|| integral gives a factor 
S{po — ko) and the vertex correction is 0 
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where u = {po+qo)/2. The logarithmic divergence was re¬ 
moved by adding the contribution from the four-fermion 
vertex in the zero sound channel. If the coupling is weak 
the scale inside the logarithm is again determined by elec¬ 
tric gluon exchange. We find that in this case the scale 
is equal to the scale in the quark self energy. 

The cancellation that occurs in the one-loop diagram 
repeats itself at any loop order if additional gluon ladders 
are added. This implies that ladder diagrams have to be 
summed. We also note that quark propagators in the lad¬ 
ders are sensitive to the small scale uj and thereby the full 
fermion self energy has to be included. A detailed anal¬ 
ysis shows that all other corrections like crossed ladders, 
vertex corrections, interconnections between the gluon 
ladders, etc. introduce extra transverse momenta and 
follow the scaling relations in the magnetic regime El 

The same phenomenon occurs in the quark-quark scat¬ 
tering amplitude. Consider a one-loop correction to the 
scattering amplitude in the BCS-channel, see Fig.|21 


^^BCsiP^g) = g^Cr jj 


^rBCs(P>^) J {k + p) S-ff {—k + p), (12) 

Here piq are the momenta of the outgoing fermions and 
Cr is a color factor. In the color anti-symmetric channel 
Cr = —{ Nc + l ) / { 2 Nc ). Combining the quark propagators 
gives 


Sij{k+p) S_i;{-k+p) 


gg (k+p) - g.g (-fc -I- p) 

2ko + Sz\p)-SzUp)-"-^f^ 


This denominator becomes sensitive to the energy ko 
when the external quarks are on-shell and the transverse 
momentum of the pair vanishes. We get 


nrlcsip. ?) - tIjJ X (s^) ■ 

(13) 

which again depends to leading order on uj = pq. The 
cutoff dependence was removed by adding the one-loop 
graph with the BCS four-fermion operator. We can com¬ 
pute the scale in the weak coupling limit. In the quark 
self energy and vertex correction the logarithmic UV di¬ 
vergence cancels between electric and magnetic gluon ex¬ 
changes and the scale is determined by transverse mo¬ 
menta k± ~ m. In the BCS channel the two terms add 
and the logarithm is sensitive to larger scales k± ^ p. 
Matching to the full tree level scattering amplitude at 
this scale gives ABCS = 2^^/^/r®/(7rTO®) as shown in da 
The double logarithmic structure of the quark-quark 
scattering amplitude in QCD was discovered by Son j^. 
Again we find that all planar, abelian ladder diagrams 
contribute at the same order but other corrections are 
perturbative and follow the scaling rules in the magnetic 
regime. The ladder sum can be determined by solving 
a Bethe-Salpeter equation, see Fig. 13 Since the quark 
propagators are sensitive to the low energy scale the 
quark self energy has to be included to all orders. We 
saw, however, that the one-loop result for the quark self 
energy is exact in the low energy limit. The situation 
simplifies further if the coupling is small. In this case 
we find that the low energy scale in the BCS channel is 
^ A ^ fiexp(—l/g). This means that the quark self 
energy correction never becomes large and can be treated 
perturbatively, too. 
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FIG. 3: The complete dynamics contributing to the Dyson- 
Schwinger equation for the BCS kernel in the low energy limit. 
The dark blob represents the full scattering amplitude and the 
shaded blob the full propagator. An analogous but indepen¬ 
dent equation is obtained for the forward channel. 

The analysis in the forward scattering channel is iden¬ 
tical to the case of an external current. Green functions 
with more than four external quark lines do not have 
kinematic enhancements except for those that occur in 
four-particle reducible graphs. The reason is that when¬ 
ever a quark line is connected to a given graph via a gluon 
an extra transverse momentum is introduced which leads 
to suppression factors at low energy consistent with the 
scaling rules. 

In summary, we have shown that QCD Green functions 
in high density QGD have a systematic low energy expan¬ 
sion in powers of The results are valid in any 

phase in which the transverse components of the gauge 
field are dynamically damped, Imn(A:o,fc) ~ m^fco/I^I- 
The damping scale m determines the asymptotic behav¬ 
ior of all gluonic correlation functions and can be used to 
define m non-perturbatively. We note, in particular, that 


the low energy interaction between gluons is weak and 
no magnetic screening mass is generated in the normal 
phase. Our results apply likewise to ordinary condensed 
matter systems described by a strong effective gauge in¬ 
teraction. 0. 

The low energy expansion is directly applicable to the 
normal phase of dense quark matter, and to modes that 
remain ungapped below the (largest) critical tempera¬ 
ture. Such ungapped modes play a vital role for a po¬ 
tential detection of quark matter in compact stars. The 
effective theory can be used to study transport proper¬ 
ties, neutrino emissivities 0 , etc. Under conditions ap¬ 
propriate to neutron star cores we have (T//r)^^^ ~ 1/10 
and the expansion is expected to converge well. This 
provides robust QCD results on the generic temperature 
dependence of low energy observables which due to the 
huge ratio /r/T are not impaired by the ignorance of low 
energy constants of order one. 

For T < Tc ~ A dense quark matter becomes super¬ 
fluid, and the Fermi liquid description breaks down com¬ 
pletely. For energies less than the gap, the non-Fermi 
liquid EFT has to be matched against an effective theory 
of the superfluid phase |l 9|| . The approximation scheme 
described here can be used to compute the parameters of 
the superfLuid EFT in an effective expansion controlled 
by parameters of magnitude (A/m)^/^ and 1/ log(m/A). 
We note, in particular, that g decreases only logarith¬ 
mically with g, while the low energy expansion contains 
inverse powers of fi. Therefore, this expansion represents 
a significant improvement over the mere use of pertur¬ 
bation theory in g and extends the range of validity of 
perturbative QCD to considerably lower densities. 
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